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Abstract 


Heating  of  continuous  carbon-fiber-reinforced  polymers  (CFRP’s)  by  the  application  of 
an  alternating  magnetic  field  has  been  shown  to  be  due  to  dielectric  losses  in  the  polymer. 
Models  that  predict  thermal  generation  in  these  composites  are  input  to  a  finite  element 
heat-transfer  analysis,  providing  the  predicted  transient  thermal  profile  in  the  plane  of  the 
laminate.  The  validity  of  the  global  thermal  generation  model  is  established  through  an 
experimental  test  matrix  in  which  various  specimen  configurations  are  evaluated  and 
compared  with  theoretical  predictions  of  transient  surface  temperatures. 
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1.  Introduction 


A  new  theory  of  dielectric  mechanisms  for  induction  heating  of  continuous  carbon-fiber 
composites  was  reported  previously  [1].  Based  on  these  mechanisms,  a  global  model  was 
developed  [2]  that  systematically  defined  the  electromagnetic  response  of  cross-ply  laminated 
composites  to  alternating  magnetic  fields  with  some  component  of  flux  normal  to  the  surface  of 
the  composite  material.  The  global  model  consists  of  three  fundamentally  separate  submodels 
that  consider  the  in-plane  electromagnetic  response  [3],  the  through-thickness  response  [4]  and, 
finally,  the  global  generation  of  heat  and  the  quantification  of  the  surface  temperature  profile 
presented  in  this  paper.  A  thermal  submodel  is  required  for  the  comparison  of  predicted  and 
measured  temperatures  at  any  point  on  the  surface  of  a  sufficiently  thin  laminated  composite 
plate.  Comparison  of  experimental  with  predicted  surface  temperatures  can  establish  the  validity 
of  the  postulated  mechanisms  of  induction  heating  in  continuous  carbon-fiber  composites. 

Figure  1  shows  an  outline  of  the  thermal  submodel  with  the  input  parameters  supplied  by  the 
foregoing  submodels  [1-5].  The  planar  grid  submodel  provides  a  two-dimensional 
(2-D)  potential  difference  distribution  for  each  ply-ply  interface.  The  fiber  layer  submodel 
provides  a  key  parameter  y  describing  the  effective  capacitive  layer  interaction  distance  and, 
thus,  the  profile  of  heating  through  the  thickness  for  each  orthogonal  or  off-axis  ply-ply 
interaction  in  the  laminate.  These  models  can  be  combined  to  generate  a  three-dimensional 
(3-D)  heat-source  distribution  in  the  laminate.  In  turn,  this  heat-source  distribution  can  be  used 
as  input  into  a  finite  element  heat-transfer  model,  which  determines  the  surface  thermal  history 
of  the  laminate.  This  latter  process  is  herein  termed  the  “thermal  submodel.” 

The  predicted  transient  surface  temperatures  obtained  from  the  thermal  submodel  are 
compared  to  experimental  surface  temperature  measurements.  ICI  Fiberite’s  prepreg  product 
APC-2  was  used  in  the  experimental  studies.  This  material  consists  of  ICI’s 
poly-etheretherketone  (PEEK)  semicrystalline  thermoplastic  resin  and  continuous  AS -4  graphite 
fibers. 
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Figure  1.  Outline  of  Thermal  Submodel. 


2.  Thermal  Conduction  Considerations 


The  heating  patterns  in  the  plane,  predicted  by  the  planar  grid  submodel  [2,  3],  and  those 
through  the  thickness,  predicted  by  the  fiber  layer  submodel  [2,  4],  combine  to  form  a  complex 
3-D  distribution  of  thermal  generation  in  the  laminate.  Although  the  models  can  predict  this 
energy  generation  point  by  point,  it  is  difficult  to  translate  these  results  into  time-variant 
temperature  distributions.  Gradients  in  the  plane  can  be  accounted  for  through  a  discretization  of 
the  surface. 

In  general,  a  3-D  transient  anisotropic  ply-ply  thermal  analysis  would  be  required.  In  our 
experiments,  relatively  thin  laminates,  where  thermal  gradients  through  the  thickness  are  shown 
to  be  negligible,  are  studied  to  confirm  the  planar  distribution  of  heating.  To  establish  this  fact, 
the  through-thickness  heat-transfer  analysis  considers  the  heat-generation  profiles  that  were 
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formulated  in  earlier  work  [2,  3],  which  can  be  somewhat  complex.  Each  possible  ply-ply 
interaction  has  such  a  profile  and  the  resulting  profile  for  a  multi-ply  laminate  may  consist  of 
several  extremum  as  illustrated  in  Figure  2. 


Figure  2.  Superposition  of  Capacitive  Layer  Interaction  Parameter,  yd* ,  Distributions  for 
[0/902]  and  [902/0]  Configurations.  The  Summation  of  These  Distributions 

Provides  the  Gamma  Distribution  for  the  [0/90]s  Configuration. 

The  profile  within  each  ply  can  be  represented  by  a  third-order  polynomial.  Steady-state  and 
transient  solutions  are  obtained  readily  for  a  single  ply-ply  interaction  using  the  equations  of 
heating,  q(x),  as  a  function  of  through-thickness  dimension  and  a  surface  convection  boundary 
condition.  If  the  polynomial  is  of  the  form  a  +  bx  +  cx2  so  that, 
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q(x)  =  P(a  +  bx+  ex?), 


(1) 


with 


0  = 


0)eoKtanSV2 


(2) 


where  to  is  the  angular  frequency,  e0  is  the  permittivity  of  vacuum  (8.85  x  10" 12  f/m),  k  is  the 
relative  dielectric  constant  of  the  polymer,  tan <5  is  the  imaginary  part  of  the  complex  dielectric 
constant  of  the  polymer,  V  is  the  fiber-fiber  interaction  potential  difference,  and  L  is  the 
characteristic  dimension  through  which  the  electric  field  created  by  V  acts  [1].  The  transient 
surface  temperature  can  be  written  as 


T(x,t)  =TX  +  yr(x,t)  +  <J>(x). 


(3) 


The  term  ¥(x,r)  is  the  transient  part  of  the  solution,  and  O(jc)  is  the  steady-state  part,  viz., 
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where  L  is  the  half-thickness  of  the  two-ply  laminate,  k33  is  the  thermal  conductivity  through  the 
thickness,  h  is  the  combined  convective  and  radiative  heat-transfer  coefficient,  and  *  is  the 
distance  from  the  ply-ply  interface  toward  the  surface  of  the  laminate.  The  transient  part  of  the 
solution  is  written  as 


V(x,t)=Y.K 
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The  steady-state  solution  can  be  obtained  independently  in  terms  of  the  temperature  at  the 
surface,  Ts: 


T=Ts+^-(L-x) 


—  [  fourth  order  polynomial  in  x] , 
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(7) 


where 

Ts  =  equilibrium  temperature  at  the  surface,  and 
Q  =  TS-  7„. 

Using  the  previous  equations  and  the  parameters  of  Table  1,  the  maximum  temperature 
difference  in  a  representative  laminate  can  be  determined  from  both  the  transient  and  steady-state 
solutions.  Figure  3  displays  the  results  for  the  representative  laminate  at  various  time 
increments.  For  the  transient  solution,  the  largest  difference  in  temperature  between  the  ply-ply 
interface,  T0,  and  the  laminate  surface,  Ts,  occurs  at  small  time.  At  time  t  =  1  s.,  this  difference  is 
approximately  5.6%  of  the  surface  temperature,  Ts.  At  t  =  «*»,  the  transient  solution  is  equivalent 
to  the  steady-state  solution  and  the  inside- to-outside  difference  is  only  0.15%  of  the  wall 
temperature. 

Due  to  the  time  increments  used  in  the  experimental  testing,  5-10  s,  and  the  accuracy 
limitations  of  the  equipment,  temperature  differences  through  the  thickness  are  not  measurable. 
Several  tests  were  performed  with  thermocouples  embedded  at  various  points  through  the 
thickness  of  a  [03/125-pm  film/903]  laminate.  No  differences  in  temperature  through  the 
thickness  were  observed.  In  contrast,  large  variations  in  the  plane  are  predicted  enabling  the 
verification  of  the  planar  grid  model  and  the  average  heating  through  the  thickness,  \q(x). 
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Table  1.  Input  Parameters  for  the  Representative  Laminate  Case  Study.  Some  Values  Are 
“Rational  Estimates”  for  Determining  the  Maximum  Possible  Temperature 
Difference  Within  a  Six-Ply  Laminate,  With  an  Interply  Resin  Thickness  of 
h0  =  127  pm 


Parameter 

Symbol 

Value 

Dimension 

Curve  Fit  Parameter 

a 

0.5169 

1/m 

Curve  Fit  Parameter 

b 

-1821.4 

I/n? 

Curve  Fit  Parameter 

c 

1.484e6 

See  Equation  (2) 

p 

23.32 

Angular  Frequency 

CO 

2n  (70,000) 

1/s 

Emissivity  Constant 

£o 

8.85e- 12 

f/m 

Dielectric  Constant 

K 

3.6 

dim.less 

Loss  Tangent 

tan5 

0.16 

dim.  less  j 

Maximum  Voltage 

V max 

68.0 

V 

Convection  Coefficient 

h 

9.1 

■■E2ZS33SHM 

Thermal  Conductivity 

kn 

0.3 

W/m°C 

Ply  Half-Thickness 

L 

4.445e-4 

m 

Coil  Current 

I 

30 

amp 

Ambient  Temperature 

25 

°c 

These  results  are  also  predicted  by  the  Biot  number,  hL/k,  and  the  rule-of-thumb  that  if 
Bi  <  0.1,  then  a  lumped-parameter  analysis  can  be  used  since  no  significant  temperature 
gradients  exist  through  the  material  of  characteristic  dimension  L.  In  this  example  for  a  six-ply 
laminate,  the  Biot  number  is  approximately  0.012.  Correspondingly,  an  APC-2  composite 
laminate  would  need  to  be  greater  than  about  16  plies  thick  (Le.,  approximately  0.2  cm  or 
0.08  in)  before  a  more  complex  3-D  analysis  would  be  required  and  appreciable  differences  in 
through-thickness  temperature  could  be  experimentally  measured.  The  following  section 
discusses  how  the  output  from  the  previous  submodels  is  consolidated  and  used  in  a  finite 
element  program  to  model  transient  surface  temperatures. 
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Figure  3.  Results  of  Closed-Form  1-D  Heat  Transfer  Solution  for  an  Example  Six-Ply 
Laminate.  T0  is  the  Interface  Temperature,  and  Ts  is  the  Surface  Temperature. 
Note  That  the  Difference  in  Transient  Temperature  Between  the  Surface  of  the 
Specimen  and  the  Center  Is  Negligible. 


3.  Elemental  Heat  Generation  and  Finite  Element 
Heat-Transfer  Solution 


The  desired  input  for  the  thermal  model  is  the  heat  generation,  qy,  in  the  planar  grid  element 
(i,j)  and  its  combined  convective  and  radiative  losses,  hy.  The  nondimensional  output  of  the 
planar  grid  submodel’s  program  for  node  “r,j”  is  A„,  where 
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V, 


(8) 


rs  r  * 

0)<pB 

The  quantity  vrs  is  the  nodal  voltage  in  volts  at  node  (r,s),  o)  is  the  angular  frequency,  and  $b 
is  the  magnetic  flux  in  webers.  The  volume-normalized  nondimensional  output  of  the  fiber  layer 
submodel’s  program  is  T: 


r  = 


y  d 

total 

b  +  n  +  m 


(9) 


where  7 total  is  the  effective  parameter  of  heating  through  the  thickness  of  the  laminate  defined  by 
the  capacitive  layer  interaction  calculations,  ct  is  the  effective  fiber  layer  separation  distance,  b  is 
the  ratio  of  the  ply-ply  separation  distance  h0  to  the  effective  fiber  layer  separation  distance  dt 
(i.e.,  b  =  hjd*),  and  m  and  n  are  the  number  of  fiber  layers  above  and  below  the  ply-ply  interface 
respectively.  The  nondimensional  parameters  Ars  and  T  are  used  to  obtain  the  total  heating  per 
element  in  the  thermal  model’s  grid.  The  grid  used  in  the  thermal  model  is  equivalent  to  that 
used  in  the  planar  grid  model.  The  nodal  voltage  in  equation  (8),  vrs,  can  be  solved  for  in  terms 
of  the  program  output,  Ars.  Since  the  thermal  model  requires  elemental  heating  values,  the  nodal 
voltages  are  averaged  over  the  element  they  represent: 


vij  =7(^.  +  ^+i+H+u  +  *Wi)» 


(10) 


where  Vy  is  the  elemental  voltage  for  the  element  (i,J)  of  the  models’  grid  (the  commas  are  used 
only  to  separate  the  subscripts).  The  elemental  heat  generation  in  watts/cubic  meter  can  now  be 
written  as, 


%  =  coeoKtan5TV^ 


(11) 
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where  the  imaginary  part  of  the  complex  dielectric  constant,  tan <5,  is  a  function  of  temperature 
[2].  From  these  procedures,  the  quantity  Xy  can  be  obtained: 


% 

tan  <5  ’ 


(12) 


with  Xy  in  watts/cubic  meter. 

The  finite  element  analysis  (FEA)  program  FEHT  was  used  for  the  heat-transfer  analysis.  A 
7x7  element  grid  was  used  as  illustrated  in  Figure  4.  A  grid  size  of  7  x  7  is  accurate  for 
centered-coil  tests  in  which  the  coil’s  flux  area  covers  at  least  14%  of  the  surface  of  the  specimen 
as  determined  by  a  convergence  study  for  the  planar  grid  model  [2],  with  an  error  of  about  7% 
from  the  infinitely  fine  grid.  This  grid  size  was  the  maximum  size  that  could  be  used  for  the 
experimental  arrangement  described  in  the  following  section.  It  represents  a  practical  choice  in 
terms  of  computer  time  limitations,  geometric  feasibility,  and  sufficient  surface  temperature 
measurement  accuracy  for  comparison  to  the  predicted  transient  surface  temperature  profile. 
Each  of  the  49  planar  grid  model  elements  from  (i,j)  =  (1,1)  to  (7,7)  were  broken  up  into  eight 
FEA  elements,  for  a  total  of  392  FEA  elements,  and  225  FEA  nodes. 

Figure  5  shows  how  each  square  model  element  was  partitioned  into  its  FEA  elements.  The 
quantity  Xy  from  equation  (12)  represents  the  thermal  generation  in  the  model  element  (ij).  Xij 
has  a  gradient  over  the  model  element  area.  This  gradient  can  be  calculated  from  the  original  vrs 
quantities  output  from  the  ply-ply  interaction  model  as  defined  in  equation  (8).  The  gradients 
between  the  four  nodal  vrs  quantities  for  each  model  element  are  assumed  to  be  linear  so  that  the 
heating  for  each  FEA  element  can  be  calculated  through  a  simple  linear  regression  involving  the 
nodal  voltage  parameters  vrs  and  X,y.  Figure  5  defines  the  nodal  locations  “a”  through  “d”  in 

addition  to  the  vrs  nodes. 
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Figure  4.  Grid  Used  in  Final  Experimental  Test  Matrix.  The  Planar  Grid  Element 
Labeled  1-10  Represent  the  Entire  Surface  Due  to  Symmetry  in  the  Cross-Ply 
Specimen  and  in  the  Central  Placement  of  the  Coil. 


Figure  5.  Mapping  From  Voltage  Surface  to  Finite  Element  Model  Nodes. 
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Assuming  linear  gradients  along  the  edges,  the  nodal  heating  values  at  all  nodes  can  be 
calculated  as 


y  V- 

y 


v  +v'~-v". 

T*  2 


b..=^ 

'  V, 


VrS+ 1  + 


Vr+ls+l  ^r^+1 


e-3 

'  v. 


..  .  ^r+U  ^  J 

”  2 


,  and 


y  V. 


y  l. 


II  .  ^r+lj+1  ^V+lj 

vr+u+  2 


(13) 


The  eight  FEA  element  heating  values  per  grid  element  (i,j)  are  given  by 


Xw  =' 


X, 

—  V^+dtj+Cy 

V  ij 


Xy,2  —  2  [Xy  +  aij  +  cij  ]* 


X,3  =- 
y.3  3 


x„ 

—Vr,s+l  +  Xy  +«j/ 


X  -I 

X,7.4  -  3 


X,. 

—vrs  +X;,  +&, 

y  y  y 


x..  =I 

y’5  3 


X, 

TTVr+U  +  Xy  +  Cy 

,  y 
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xii6  =  - 

y,6  3 


X„ 

T^'Ur+l,s+Xij  +  dij 

.  y 


Xij,7  =  ~  [xy +  by  +  dy  ],  and 


X, 

vsl 


y+i.s+i +bij+dij 


(14) 


The  average  of  Xy,;  through  X,y,§  is  not  exactly  equal  to  Xy  due  to  the  assumption  of  linear 
gradients  between  adjacent  nodes;  however,  this  error  was  consistently  less  than  the  7%  error 
inherent  in  the  planar  grid  analysis  for  the  specimens  modeled  and  tested.  For  the  centered-coil 
tests  used,  some  model  elements,  (i,j),  had  equivalent  heat-generation  input  values.  Due  to  this 
symmetry,  10  distinct  input  values  were  possible,  as  shown  in  Figure  4.  The  “heat  generation” 
term  needed  for  the  thermal  calculation  is 


qu  =  tan  8(t)~ 


2  hu{T) 


(15) 


where 

<lij,k  -  total  heat  flux  in  FEA  element  k  of  model  element  (ij), 

Xij'k  =  output  after  linear  regression  fits, 

tan 5(7)  =  dielectric  loss  tangent  as  a  function  of  element  temperature, 

hj,k  =  combined  convective  and  radiative  heat  transfer  losses  from  each  surface  of  the  plate  as  a 
function  of  element  temperature,  and 

L-(b+m  +  n)a  =  thickness  of  plate. 
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The  two  fundamental  parameters,  tan<5(7)  and  %*(7),  were  experimentally  investigated.  A 
study  of  the  loss  tangent  properties  of  PEEK  was  performed  and  is  presented,  in  detail,  in 
Fink  [2]  and  Fink,  McCullough,  and  Gillespie  [5].  The  convection  heat  transfer  coefficient,  h, 
was  experimentally  determined  for  the  test  configuration.  This  analysis  is  discussed  in  Fink  [2]. 

Other  input  to  the  program  includes  the  density,  heat  capacity  and  thermal  conductivity  of  the 
material.  The  FEA  program  requires  the  input  of  effective  quasi-isotropic  thermal  properties. 
The  thermal  conductivity  of  an  APC-2  prepreg  ply  is  highly  anisotropic;  the  value  in  the  fiber 
direction,  is  approximately  6.20  W/m°C,  and,  in  the  transverse  direction,  kj-i,  is 
approximately  0.34  W/m°C  from  fiber  volume  fraction  calculations  using  polymeric  and 
anisotropic  fiber  properties  (Table  2)  and  a  fiber  volume  fraction  of  0.61.  For  orthotropic  cross- 
ply  laminates  in  which  the  number  of  0°  plies  is  equal  to  the  number  of  90°  plies  (M  =  N), 
isotropic  values  for  thermal  conductivity  can  be  used. 

Table  2.  Material  Properties  for  Poly-Etheretherketone  and  AS-4  Carbon  Fibers 


Specific  Heat 
J/kg°C 

Thermal  Conductivity 
W/m°C 

Matrix  (PEEK)  [6] 

1,340 

1,320 

l—  0.2512 

Fiber  (AS-4  Graphite)  [8] 

921 

1,800 

0.426  (trans.) 

10.0  (long.) 

Consider  a  [0A///io_thick_fihn/90w]  cross-ply  laminate  consisting  of  M  0°  plies,  an  interlayer 
resin  thickness  of  h0,  and  N  90°  plies.  A  “longitudinal”  value  for  the  thermal  conductivity,  L,  of 
the  laminate  can  be  estimated  through  volume  averaging: 


where 


j^la  miniate 


ML?  +  Nk?  +  -^k 

n  22  127  , 


22 


M  +  N  + 


127 


(16) 
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km  =  thermal  conductivity  of  matrix  «  0.2512  W/m°C,  and 


h0  is  measured  in  microns. 

For  M  -  N  =  \  and  h0  -  8  |Ltm  (a  normally  consolidated  [0/90]  cross-ply  laminate),  ki  for  the 
laminate  is  approximately  3.18  W/m°C.  For  balanced  laminates,  M  =  N,  the  longitudinal  and 
transverse  conductivities  are  equivalent,  ki  =  kj.  As  long  as  thermal  gradients  through  the 
thickness  are  minimal,  the  in-plane  thermal  conductivity  is  isotropic.  Similar  calculations  are 
possible  for  determining  the  density,  p,  and  specific  heat,  cp,  values  for  the  various 
configurations  under  consideration. 

4.  Experimental  Procedures 


The  thermal  submodel  described  previously  completes  the  global  model  of  heating  in 
cross-ply,  continuous  carbon-fiber-reinforced  composites.  The  local,  in-plane,  and 
through-thickness  submodels  are  combined  to  provide  input  to  a  finite  element  heat-transfer 
analysis,  providing  the  predicted  transient  thermal  profile  in  the  plane  of  thin  laminates.  Error 
calculations  concerning  through-thickness  thermal  gradients  and  linearity  of  potential  field 
gradients  in  the  plane  were  considered  in  the  solution.  As  shown  later,  the  validity  of  the 
assumptions  that  encompass  the  thermal  submodel  are  supported  by  the  close  agreement  with  the 
model’s  output  with  experimental  data. 

The  following  sections  provide  experimental  evidence  supporting  the  global  model.  The 
various  submodels  previously  discussed  [1-5]  established  the  significant  physical  parameters 
required  to  characterize  the  response  of  cross-ply  composites  to  alternating  magnetic  fields. 
Various  material  and  experiment-related  quantities  are  evaluated  individually  and  tabulated  with 
their  corresponding  errors.  An  experimental  test  matrix  was  performed  that  represents  a  practical 
range  of  the  significant  parameters  and  the  effects  of  changing  those  parameters  individually. 
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Table  3  shows  the  20.3-cm-square  (8  in)  plate  lay-up  configurations  fabricated  and  tested 
using  coil  currents  of  6,  10,  and  20  amp;  the  matrix  consists  of  75  separate  tests.  The  placement 
of  the  10  surface  thermocouples  is  shown  in  Figure  4.  Specimens  were  labeled  as  indicated  in 
Table  3  and  are  referred  to  by  the  codes  shown  in  parentheses,  e.g.,  (0_11)  represents  the  0-mil 
interlayer  [0i/90i]  specimen  and  (5_22)  represents  the  5  mil  interlayer  [O2/9O2]  laminate 

specimen  where  1  mil  =  25.4  pm.  Specimens  are  further  referred  to  as,  for  example,  (0_1 1 _ 6), 

where  the  “6”  represents  the  coil  current  used  in  amps. 


Table  3.  Laminate  Configurations  Processed  for  Final  Test  Matrix.  Each  20.3-cm-Square 
(8  in)  Specimen  Was  Tested  With  a  10.2-cm-Diameter  (4  in)  Centered  Helmholtz- 
Type  Coil  at  6-,  10-,  and  20-amp  Applied  Coil  Current. 


Nominal  Interlayer  PIIK 

Resin  Film  Added,  h0  (1  mil  =  25.4  pm)  | 

(None) 

1  mil 

2  mil 

3  mil 

5  mil 

(1,1) 

A 

(0_H) 

B 

(1_11) 

C 

(2_11) 

D 

(3_11) 

E 

(5_11) 

(2,2) 

F 

(0_22) 

G 

(1_22) 

H 

(2_22) 

I 

(3_22) 

J 

(5_22) 

(M,N) 

(3,3) 

K 

(0_33) 

L 

(1_33) 

M 

(2_33) 

N 

(3_33) 

O 

(5_33) 

(1,2) 

P 

(0_12) 

Q 

G_12) 

R 

(2_12) 

S 

(3_12) 

T 

(5_12) 

(1,3) 

U 

(0_13) 

V 

(1_13) 

W 

(2_13) 

— 

Y 

(5_13) 

DuPont  Co.  provided  the  use  of  an  Ameritherm  50-200-kHz  15-kW  induction  generator. 
The  machine  self-tunes  the  frequency  for  the  coil  and  load  placed  on  it.  Output  includes  the 
frequency,  current  to  coil,  and  voltage  applied.  The  current  to  the  coil  reading  has  about  a  2-amp 
error  in  the  reading  itself,  which  can  be  100%  of  the  total  reading  when  applying  low  power  to 
the  coil.  The  induction  coils  were  fabricated  using  1/4-in-outside-diameter.  Copper  tubing  and 
compression  fittings  using  a  Helmholtz  design,  which  provides  a  very  uniform  field  with  all  field 
lines  normal  to  the  plane  of  the  specimen  and  with  constant  flux  across  the  plane  of  the 
specimen.  The  coil  design  used  in  our  experiments  tunes  to  the  Ameritherm  equipment  at 
70  kHz.  The  20.3-cm-square  (8  in)  specimens  were  placed  in  a  wooden  frame,  which  serves  two 
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purposes.  First,  it  grips  the  specimen  along  all  four  edges  (with  1/8-in  overlap  at  the  edges)  and 
keeps  the  unbalanced  (1,2)  and  (1,3)  and  the  unsymmetrical  laminates  (all  specimens)  square  in 
the  plane.  Second,  the  frame  assembly  slides  vertically  between  two  grooved  wooden  guides 
attached  to  a  wooden  base,  allowing  the  frame  and,  thus,  the  specimen  to  be  rapidly  dropped  into 
the  stationary  Helmholtz  coil  position.  This  provides  an  instant-on  capability  for  the  field  not 
otherwise  possible  with  the  Ameritherm  induction  equipment  at  such  low  coil-currents,  due  to 
the  self-tuning  nature  of  the  equipment.  The  thermocouples  were  placed  on  the  specimen  and 
each  read  by  a  digital  E-type  thermocouple  reader  (a  computerized  input  was  not  available).  The 
accuracy  of  these  readers  is  ±1  °C.  The  process,  including  the  thermocouple  readers  and  a  digital 
timer,  was  recorded  on  videotape,  which  was  later  reviewed  to  record  the  thermocouple  readings 
at  10-s  increments.  These  time  increments  provided  sufficient  accuracy  for  comparison  to 
predictions. 

5.  Model  Predictions 

Table  4  indicates  the  values  used  in  the  models  for  comparison  with  the  final  experimental 
matrix.  The  application  of  the  predictive  models  is  illustrated  for  the  laminate  configuration 
(0_1 1_6).  Nondimensional  laminate  and  induction  source  parameters  (generally,  those  without 
error  values)  were  input  into  the  planar  grid  in-plane  interaction  submodel,  providing  a 
dimensionless  nodal  voltage  matrix.  Separately,  dimensional  laminate  characteristics  and  their 
respective  errors  were  input  to  the  fiber  layer  through-thickness  interaction  submodel,  providing 
the  total  effective  inverse  fiber-layer  interaction  distance  per  unit  thickness.  These  outputs  were 
combined  with  additional  material  and  laminate  properties  for  input  into  a  program  based  upon 
the  local  and  in-plane  interaction  submodels.  This  output  provides  the  7  x  7  heat-generation 
matrix  per  unit  loss  tangent  with  upper  and  lower  error  bounds.  These  results  were  then  used  as 
input  to  the  2-D  FEA  heat-transfer  submodel  along  with  loss  tangent  and  surface  heat-loss 
functions  with  temperature  to  provide  the  predicted  transient  thermal  behavior  at  the  surface  of 
the  laminates. 
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Table  4.  Values  and  Associated  Errors  for  Material  and  Experimental  Parameters 


Property 

Symbol 

Source 

Coil  Current 

/ 

Machine 

Reading 

Magnetic  Flux 

4>b 

Coil  Calc’s 
[2] 

Frequency 

/ 

Machine 

Reading 

Dielectric  Constant 

e 

Measured 

[5] 

Dielectric  Loss  Tangent 

tan  S 

Measured 

[5] 

Fiber  Diameter 

df 

[8] 

Fiber  Volume  Fraction 

xf 

[9] 

Radiation  Heat-Transfer 
Losses 

R 

Estimated 

[2] 

Number  of  Fiber  Layers 

Per  Ply 

m/M,n/N 

Measured 

[2] 

lActual  Interply  Thickness 

ho/Nom. 

Measured 

fPer  Unit  Nominal  Input 

Input 

[2] 

Total  Laminate  Thickness 

2L 

Calculated 

Wall  Temperature 

Ts 

Measured 

Convective  Heat-Transfer 
Coefficient 

h 

Measured 

[2] 

Heat  Capacity 

Cp 

Calculated 

Density 

P 

Calculated 

Thermal  Conductivity 

k 

Calculated 

Value 

Error 

Units-  j 

6,  10, 20 

±2 

Amps 

1.219e-6*I 

Lo:  1.08e-6*I 
Hi:  1.2196-6*1 

Webers 

70,000 

±2000 

Hertz 

3.72 

±0.10 

dim.less 

Varies  With  T 
[2] 

3%  ' 

[2]  I 

dim.  less 

8 

±0.5 

microns 

0.61 

±0.01 

dim.less 

Varies  With  T 
[2] 

N/A 

17.3 

±  1.6 

dim.less 

HI 

dim.less 

Varies  With 
Specimen 

Varies  With 
Specimen 

meters 

Varies  (Depends 
on  Rate) 

±1 

■l 

9.1 

±0.1 

[2] 

N/A 

J/k°C 

[2] 

N/A 

m 

[2] 

N/A 

W/m°C 

Equation  (17)  is  the  input  unit  magnetic  flux  matrix  for  the  planar  grid  model  (the  total  input 
sums  to  unity).  This  input  best  represents  the  shape  and  placement  of  the  coil  and  the 
distribution  of  flux  on  the  surface  for  a  7  x  7  grid.  Note  that,  although  the  distribution  of  flux  is 
continuous  in  the  actual  laminate  plane,  the  coarse  modeled  grid  of  the  planar  grid  model 
necessarily  has  some  partially  filled  elements  as  has  been  previously  established  [2]. 


».] 

it*.. 


j= 1  <=1 


0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0.071 

0.143 

0.071 

0 

0 

0 

0 

0.143 

0.143 

0.143 

0 

0 

0 

0 

0.071 

0.143 

0.071 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

(17) 


Equation  (18)  displays  an  example  of  the  8  x  8  matrix  of  dimensionless  nodal  voltage  (i.e., 
voltage  per  unit  frequency  per  unit  flux)  for  the  (0_11_6)  laminate  analysis.  The 
nondimensionalized  total  effective  inverse  interaction  distances,  y,  for  the  predicted  parameters 
and  their  errors  (Table  4)  for  each  specimen  type  (and  10-mil  [254  pm]  interlayer  thickness) 
were  determined  previously  [2]. 


0:757... 

0.840 

0.943 

0.426 

0.426 

0.943 

0.840 

0.757" 

0.840 

6:943.....  1.08 

0.490 

0.490 

1.08 

0.943 

0.840 

0.943 

1.08 

OO 

!N 

0.581 

1.27 

1.08 

0.943 

0.426 

0.490 

0.081 

6:i7-5.. 

0.175 

0.081 

0.490 

0.426 

0.426 

0.490 

0.081 

0.175 

0.175 

0.081 

0.490 

0.426 

0.943 

1.08 

1.27 

0.581 

0.581 

1.27 

1.08 

0.943 

0.840 

0.943 

1.08 

0.490 

0.490 

1.08 

0.943 

0.840 

0.757 

0.840 

0.943 

0.426 

0.426 

0.943 

0.840 

0.757 

Note  in  equation  (18)  that  the  dotted  lines  delineate  symmetry  about  the  one-eighth  section. 
In  the  equations  that  follow,  this  symmetry  is  noted. 
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Several  observations  derived  from  parametric  studies  are  useful  in  predicting  the  trends  of 
through-thickness  heating  as  a  function  of  the  number  of  plies  (Mfl)  on  either  side  of  the 
interface  and  of  the  thickness  of  the  matrix-rich  interface,  h0.  Figure  6  shows  the  effect  of 
increasing  M  for  various  nominal  values  of  h0  (data  is  shown  without  error  bars  for  ease  of 
reading).  Note  that  the  effect  (i.e.,  the  change  in  heating  with  increasing  M)  is  greater  for  larger 
interface  thicknesses  due  to  the  dominance  of  the  parameter  h0.  Figure  7  shows  the  effect  of 
increasing  h0  for  the  various  laminate  configurations.  Note  that  the  (1,2)  and  (1,3)  configurations 
do  not  differ  significantly  in  heating  and  that  the  (1,3)  and  (2,2)— each  with  equivalent  numbers 
of  plies — heat  significantly  differently;  both  of  these  observations  are  due  to  the  ineffectiveness 
of  increasing  the  number  of  plies  on  one  side  of  the  interface.  Finally,  it  is  of  interest  to  compare 
the  predicted  overall  “heating”  of  the  fiber  layer  interaction  model  with  a  situation  in  which 
dielectric  heating  is  restricted  to  the  interface  region,  h0.  Figure  8  compares  the 
nondimensionalized  “effective”  inverse  interaction  distances  from  the  present  model  with  the 
nondimensionalized  inverse  ha  (note  the  inverse  x-axis  scale  to  represent  increasing  interface 
thickness).  As  the  interface  thickness  gets  larger  (i.e.,  inverse  h0  approaches  zero),  the  difference 
between  the  two  models  approaches  zero  since  h0  dominates  the  heating  at  low  values. 
However,  at  practical  interface  thicknesses  (even  for  specimens  with  no  additional  PEEK 
interlayer  film),  the  present  model  predicts  heating  that  is  2-3  orders  of  magnitude  higher  than 
the  situation  in  which  heating  occurs  only  at  the  interface. 

The  planar  voltage  profile  of  equation  (18)  and  the  total  effective  through-thickness 
interaction  parameter,  y,  were  used  to  predict  the  heating  profile  in  the  plane  of  the  specimen  in 
watts/cubic  meter  per  unit  loss  tangent  as  shown  in  equations  (19a-c).  “B”  represents  the 
predicted  values,  while  “U”  and  “L”  represent  the  upper  and  lower  bounds  calculated  from  the 
values  and  errors  tabulated  in  Table  4. 
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2.5 


2 


gamma  5 
per  unit 
thickness 
(l/um2) 


0.5 


0  -I - 1 - 1 - 1 - 1 - 1— * 

0  1  2  3  4  5 

Number  of  Piles  on  Each  Side  of  Interface  (M  =  N) 

Figure  6.  Predictions  of  Capacitive  Layer  Model  Study  of  Test  Matrix  Specimens  Showing 
the  Effect  of  Increasing  Ply  Thickness  on  Heating. 
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Figure  7.  Predictions  of  Capacitive  Layer  Model  Study  of  Test  Matrix  Specimens  Showing 
the  Effect  of  Increasing  Interface  Thickness  on  Heating. 


20 


12QG0 
10QG0 
8Q00 

(gamma)d*  JQQQ 

4Q00 

2aoo 

QOO 

Q6  Q5  Q4  Q3  Q2  Q1  0 
(1 Jh0)d* 


Figure  8.  Predictions  of  Capacitive  Layer  Model  Study  of  Test  Matrix  Specimens 
Comparing  Nondimensional  “Effective”  Inverse  Interaction  Distances  With 
the  Nondimensionalized  h0.  Note  the  Inverse  x-Axis  Scale  to  Represent 
Increasing  Interface  Thickness. 
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Note  that  the  symmetrical  placement  of  the  coil  on  the  specimen  and  the  circular  shape  of  the 
coil  leads  to  only  10  unique  values  in  the  matrix  Q.  The  325%  increase  in  heating  from  the 
lower  bound  to  the  upper  bound  appears  substantial.  However,  this  difference  translates  into 
only  a  15  °C  difference  in  surface  temperature  at  the  hottest  point  on  the  surface,  as  shown  in 
Figure  9. 


Time  (seconds) 

Figure  9.  Experimental  Surface  Temperature  History  for  Thermocouple  No.  5  on  (0_11_6) 
Test  Specimen. 
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Equations  (20a-c)  show  the  predicted  surface  temperature  profile  and  the  corresponding 
upper  and  lower  bounds  respectively. 
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Large  differences  between  the  theoretical  bounds  in  surface  temperature  were  common; 
especially  for  the  lower  currents.  In  addition  to  the  large  machine-induced  error  in  current 
readings,  other  factors  add  to  the  disparity  in  the  bounds  of  the  thermal  response.  The  sensitivity 
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of  the  microstructural  properties  of  fiber  diameter  and  fiber  volume  fraction  and  the 
macrostructural  properties  of  interply  resin-rich  region  thickness  and  average  number  of  fibers 
through  the  thickness  of  a  single  ply  have  significant  effects  on  the  predicted  bounds.  However, 
experimental  results  are  later  shown  to  verify  the  predicted  solution  despite  the  wide  bounds. 
This  is  due  to  the  averaging  of  these  properties  over  millions  of  actual  nodal  locations  (fiber- 
fiber  intersections),  which  tends  to  smooth  out  the  effects  of  any  deviations  from  the  average  in 
these  structural  properties. 

6.  Comparison  of  Predicted  and  Observed  Temperature 
Profiles 

Although  experimental  results  for  the  nonorthotropic  laminates  [e.g.,  (1,2)  and  (1,3)]  were 
obtained  [2],  comparable  theoretical  results  were  not  feasible  due  to  the  FEA  program’s  inability 
to  handle  anisotropic  thermal  conductivity  properties.  Accordingly,  attention  is  directed  to  the 
orthotropic  laminates  (1,1),  (2,2),  and  (3,3).  The  model  results  are  obtained  from  the  FEA  nodal 
transient  temperatures  for  the  10  nodes  corresponding  to  the  thermocouple  locations,  as  shown 
previously  in  Figure  4.  Only  the  results  for  the  10-amp  applied  coil-current  tests  are  considered 
here.  This  current  resulted  in  significant  heating  for  all  specimens,  while  remaining  under  the 
melting  temperature  of  PEEK  and,  more  importantly,  within  the  range  of  temperatures  for  which 
the  loss  tangent  was  measured.  The  (0_1 1_6)  specimen  is  examined  in  detail.  Fink  [2]  presents 
extensive  additional  data  for  other  specimens. 

6.1  Planar  Temperature.  Figure  10  shows  the  experimental  transient  surface  temperature 
results  over  the  first  3  min.  As  noted  previously,  Figure  9  shows  the  predicted  results  for  the 
point  of  highest  heating  on  the  surface  of  the  specimen  with  the  experimental  results  for  the  point 
of  highest  heating  superimposed  with  the  appropriate  error  bars.  The  equilibrium  experimental 
value  is  45  °C,  while  the  theoretical  equilibrium  value  is  46  °C,  representing  a  (l/(46-23))xl00  = 
4%  decrease  in  temperature  gain.  Figure  11  compares  the  equilibrium  results  for  all 
thermocouples  with  the  corresponding  predicted  values.  Note  that  the  No.  2  and  No.  6 
thermocouple  positions  heated  more  than  the  theory  predicted,  while  the  No.  1  and  No.  3 
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Figure  10.  Experimental  Transient  Surface  Temperature  Results  for  Specimen  (0_11_6) 
for  All  Nodal  Positions. 


positions  heated  less  than  the  model  predicted.  Figure  12  indicates  the  thermocouple  positions 
on  the  specimen  and  the  percent  difference  of  observed  and  predicted  equilibrium  temperature 
averaged  over  all  the  10-amp  tests  and  normalized  against  the  ratio  of  experimental  to  predicted 
equilibrium  surface  temperature  for  position  No.  5.  The  standard  deviation  in  these  results  was 
less  than  2%.  Note,  again,  the  significant  increase  in  the  No.  2  and  No.  6  positions  and  the  slight 
decrease  in  the  No.  1,  No.  3,  and  No.  8  positions.  This  expected  result  is  explained  by  the  1-D 
assumption  and  the  analysis  of  effective  properties  in  which  individual  plies  are  not  modeled  in 
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Figure  11.  Experimental  and  Predicted  Equilibrium  Surface  Temperature  Data  for 
Specimen  (0_11_6)  Including  Appropriate  Error  for  Each. 

the  heat-transfer  analysis  used  in  this  study.  Position  No.  5  represents  the  point  of  highest 
heating  in  the  specimen  quadrant  shown  in  Figure  12.  The  1-D  assumption  may  break  down  in 
light  of  the  severe  thermal  gradients  around  thermocouple  position  No.  5.  An  analysis  of  the  3-D 
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Figure  12.  Diagram  Showing  Average  Variation  of  Experimental  Equilibrium  Results 
From  Predicted  Equilibrium  Results  for  All  10  Thermocouples. 

heat- source  distribution,  provided  by  the  in-plane  planar  grid  and  the  through-thickness 
capacitive  fiber  layer  models,  which  would  provide  a  layer-by-layer  anisotropic  heat-transfer 
analysis  would  likely  resolve  these  discrepancies.  This  is  a  topic  of  future  work. 

The  positions  of  Nos.  2  and  6  in  Figure  12  lie  along  fiber  paths  in  the  0°  and  90°  plies; 
therefore,  they  are  expected  to  heat  somewhat  faster  than  predicted  due  to  the  higher  actual 
thermal  conductivity  in  these  directions  vs.  the  effective  properties  used  in  our  analysis. 
Likewise,  the  Nos.  1,  3,  and  8  positions  are  expected  to  heat  slower  than  predicted.  The  higher 
values  in  the  Nos.  4  and  7  positions  are  not  meaningful  since  the  actual  temperatures  measured  in 
these  locations  was  very  low.  Position  Nos.  9  and  10  were  also  low  heating.  Their 
percent-difference  results  are  not  included  due  to  the  insignificant  temperature  values  recorded  at 
these  locations. 
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Again,  tabular  equilibrium  results  for  all  the  10- amp  applied  coil-current  tests  are  available 
[2].  These  tables,  however,  fail  to  indicate  the  transient  behavior  of  the  surface  temperature, 
which  would  provide  comparative  information  on  heating  rates.  One  way  of  providing  this 
information  is  to  represent  transient  behavior  as  a  linear  increase  in  temperature  with  time  to  an 
equilibrium  temperature,  Teq,  at  time  teq.  The  time  to  equilibrium  is  calculated  from  the  transient 
results  as  the  time  to  reach  a  temperature  within  3%  of  the  equilibrium  temperature  (taken  at  t  = 
180  s)  as  illustrated  in  Figure  13. 


Figure  13.  Schematic  Definition  of  Equilibrium  Temperature  and  Equilibrium  Time.  The 
Equilibrium  Temperature  Is  Defined  as  the  Temperature  Measured  at 
T  =  180  s.  The  Equilibrium  Time  Is  Defined  as  the  Time  at  Which  the 
Temperature  Has  Reached  3%  of  the  Equilibrium  Temperature.  The  “Time- 
Temperature  Equilibrium  Point”  Is  Marked  With  the  Small  Circle. 


Although  the  various  thermocouples  (i.e.,  surface  positions)  reached  equilibrium  at  different 
times,  the  average  time  to  equilibrium  is  valuable  for  comparison  between  specimens.  Figure  14 
shows  these  calculated  points  in  the  thermal  history  data  for  the  assumed  linear  heating  profiles. 
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Time  for  Surface  to  Reach  Equilibrium  (s) 

Figure  14.  Summary  of  Time-Temperature  Equilibrium  Points. 
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Figure  14a  shows  that,  as  expected,  Teq  and  teq  both  increase  with  increasingly  ply  thickness;  Teq 
increases  due  to  the  higher  heating  predicted  by  our  theory,  and  teq  increases  due  to  the  increased 
amount  of  material  being  heated.  Figure  14b  shows  that  increasing  the  interlayer  thickness 
significantly  decreases  Teq  as  predicted  and  slightly  increases  teq  due  to  less  significant  increases 
in  added  material;  however,  the  material  that  is  added  has  a  low  thermal  conductivity  and  is  in 
the  center  of  the  specimen.  Figure  14c  shows  that  an  increase  in  the  current  to  the  coil  increases 
Teq  dramatically,  yet  the  time  to  reach  equilibrium  does  not  change. 

6.2  Effect  of  Varying  Interface  Thickness.  The  predicted  effect  of  varying  the  interface 
thickness  on  the  parameter  y  (given  in  Figure  7)  indicated  a  rapid  decrease  in  the  heating 
parameter  y  with  increasing  interface  resin  thickness  h0.  As  noted  previously,  the  point 
corresponding  to  the  location  of  thermocouple  No.  5  in  Figure  4  represents  the  highest  predicted 
and  recorded  heating  in  the  plane.  Figure  15  shows  experimental  data  for  the  equilibrium  (steady 
state)  surface  temperature  at  this  point  for  the  [0/90],  [O2/9O2]  and  [O3/9O3]  specimens  at  10-amp 
current.  The  error  bars  on  the  experimental  results  are  indicative  of  the  ±1  °C  reading  of  the 
digital  thermocouple  reader  and  an  additional  ±5%  error  on  the  accuracy  of  the  placement  of  the 
thermocouple  on  the  specimen.  This  value  of  5%  was  determined  by  analyzing  the  changes  in 
temperature  over  one  square  centimeter  centered  on  the  point  of  highest  heating.  Figure  16 
compares  predicted  and  experimental  steady-state  surface  temperatures  at  this  point  for  the 
[0/90]  specimen  [code  (x_ll_10)].  Experimental  results  agree  quite  well  with  predictions. 
Equivalent  comparisons  were  obtained  for  other  specimen/coil  configurations  [2]. 

6.3  Effect  of  Varying  Ply  Thickness.  Figure  6  shows  the  predicted  effect  of  increasing  the 
ply  thickness  on  the  capacitive-layer-model’s  parameter  y.  Moderate  increases  in  the  heating 
parameter  y  were  predicted  for  substantial  increases  in  ply  thickness.  Figure  17  shows  the 
observed  effect  of  increasing  the  ply  thickness  on  the  equilibrium  surface  temperature  at 
thermocouple  No.  5  for  various  interply  thicknesses.  Although  increasing  ply  thickness  has  little 
effect  on  the  parameter  y,  it  has  a  more  significant  effect  on  the  temperature.  This  is  due  to  the 
significant  increase  in  loss  tangent  of  PEEK  at  the  glass  transition  temperature,  resulting  in 
increased  heating  in  those  specimens  that  approach  a  temperature  of  140  °C,  as  shown  in 
Figure  6. 
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Equilibrium  Temperature  (°C) 


Nominal  Interface  Resin  Thickness  (mil) 

Figure  15.  Experimental  Results  Showing  Effect  of  Varying  the  Thickness  of  the 
Interfacial  Resin  Layer. 
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Equilibrium  Surface  Temperature  (°C) 


□  Experiment 
■  Theory 


Nominal  Interface  thickness  (mil) 


Figure  16.  Comparison  of  Predicted  and  Observed  Surface  Temperatures  Showing  Effect 
of  Increasing  Interface  Thickness,  h0,  on  Equilibrium  Surface  Temperature. 
The  Heavy  Lines  Represent  the  Experimental  Results,  While  the  Thin  Lines 
Represent  Predicted  Results  With  Upper  and  Lower  Bounds. 


•* 
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Number  of  Piles  on  Each  Side  of  Interface  (M  =  N) 

Figure  17.  Experimental  Observation  Showing  Effect  of  Increasing  Ply  Thickness  on 
Equilibrium  Surface  Temperature  for  Various  Interface  Resin  Thicknesses. 


6.4  Effect  of  Varying  Coil  Current.  Figure  18  shows  the  effect  of  varying  the  applied  coil 
current  (i.e.,  increasing  the  magnetic  flux)  on  the  (0_11_7)  specimens  using  currents,  I,  of  6,  10 
and  20  amp.  A  current-squared  power  relationship  is  expected.  The  curve  fit  shown  on  the 
graph  is  for  the  equation  T  =  0.7(I)2  +  23  °C.  With  the  loss  tangent  changing  with  temperature, 
the  thermal  data  is  expected  to  show  a  faster  rise  in  temperature  than  for  a  constant  loss  tangent 
value.  Therefore,  the  current-squared  relationship  includes  both  the  resistive  nature  of  heating 
and  the  variation  of  loss  tangent  with  temperature,  which  is  offset  by  the  temperature-dependent 
convection  and  radiation  losses. 
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Current-to-Coil  (amp) 

Figure  18.  Comparison  of  Predicted  and  Observed  Surface  Temperatures  Showing  Effect 
of  Applied  Coil  Current  on  Heating. 


7.  Summary 


Local,  in-plane,  and  through-thickness  submodels  were  combined  to  provide  input  for  a  finite 
element  heat-transfer  analysis  to  predict  transient  thermal  profile  in  the  plane  of  thin  laminates. 
Assumptions  concerning  through-thickness  thermal  gradients  and  linearity  of  potential  field 
gradients  in  the  plane  were  imposed  to  obtain  a  solution.  The  validity  of  the  assumptions  that 
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encompass  the  thermal  submodel  are  supported  by  the  agreement  with  the  model’s  output  with 
experimental  data.  The  convection  heat-transfer  coefficient  was  experimentally  determined. 

An  experimental  test  matrix  and  apparatus  were  employed  to  verify  the  proposed  theory  of 
dielectric  heating  as  the  fundamental  mechanism  of  thermal  generation  in 
continuous  carbon-fiber  thermoplastic-matrix  cross-ply  laminated  composites  subjected  to  a 
transverse  magnetic  field.  It  has  been  shown  that,  as  the  ply-ply  interface  resin-rich  region 
thickness  h0  increases,  the  heating  in  the  laminate  decreases  significantly.  This  indicates  that 
Joule  heating  due  to  fiber-fiber  contact  is  not  likely.  Furthermore,  increasing  the  number  of  plies 
on  either  side  of  the  interface  increases  the  heating  indicating  that  thermal  generation  in  the 
laminate  is  a  function  of  ply  thickness  and  most  likely  occurs  throughout  the  thickness  of  the 
laminate. 

The  accuracy  of  the  surface  temperature  measurements  to  those  predicted  under  a  variety  of 
physical  parameters  gives  credence  to  the  validity  of  the  original  proposal  stated  in  Fink, 
McCullough,  and  Gillespie  [1].  Claims  of  Joule  losses  in  the  fibers  [10-13],  fiber  contact 
resistance  losses,  and  local-only  heating  at  the  ply-ply  interface  are  not  supported  by  the  current 
experimental  investigation.  Joule  losses  in  the  fibers  would  result  in  constant  heating  on  the 
perimeter  of  the  Helmholtz-type  coil — this  was  not  the  case;  heating  was  shown  to  occur  at 
precise  points  and  with  specific  intensities  as  predicted  by  the  global  model.  Fiber  contact 
resistance  losses,  while  potentially  exhibiting’ similar  qualitative  heating  patterns,  would  induce 
significant  current  in  the  fibers  themselves,  resulting  in  measurable  Joule  heating  in  the  fibers. 
This  was  not  witnessed  in  this  study.  Although  not  directly  shown,  the  existence  of  the  gradient 
of  heating  through  the  thickness  of  the  laminate,  as  predicted  by  the  capacitive  fiber  layer 
submodel,  is  indirectly  supported  by  the  model’s  ability  to  predict  the  total  thermal  generation. 
However,  the  quantitative  summation  of  through-thickness  heat  generation  was  verified.  If 
heating  only  occurred  at  the  resin-rich  ply-ply  interface  {h0  region),  then  the  expected  total 
heating  would  be  significantly  lower  and  result  in  significantly  lower  surface  temperatures. 
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Continuing  work  involves  extending  the  present  theory  and  models  to  laminates  containing 
off-axis  plies  and  performing  a  3-D  heat-transfer  analysis  using  anisotropic  properties  and  a  ply- 
by-ply  approach.  This  will  allow  the  investigation  of  thick  laminates  to  further  verify  the 
through-thickness  thermal  generation  model  and  provide  insight  into  special  problems  associated 
with  angle-ply  laminates. 
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D  RIGOGLIOSO 
PICATINNY  ARSENAL  NJ 
07806-5000 
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1  COMMANDER 

US  ARMY  ARDEC 

AMSTAARSRE 

DYEE 

PICATINNY  ARSENAL  NJ 
07806-5000 

1  COMMANDER 

US  ARMY  ARDEC 
AMSTA  AR  WET 
T  SACHAR 
BLDG  172 

PICATINNY  ARSENAL  NJ 
07806-5000 

1  COMMANDER 

US  ARMY  ARDEC 
SMCAR  ASF 

PICATINNY  ARSENAL  NJ 
07806-5000 

I  COMMANDER 
US  ARMY  ARDEC 
AMSTA  ARWELF 
INTELLIGENCE  SPECIALIST 
M  GUERRIERE 
PICATINNY  ARSENAL  NJ 
07806-5000 

I I  PROJECT  MANAGER 
TANK  MAIN 
ARMAMENT  SYSTEMS 
SFAE  GSSC  TMA 

R  MORRIS 
CKIMKER 
D  GUZOWICZ 
EKOPACZ 
RROESER 
R DARCY 
R  MCDANOLDS 
LDULISSE 
C  ROLLER 
J  MCGREEN 
B PATTER 

PICATINNY  ARSENAL  NJ 
07806-5000 


NO.  OF 

COPIES  ORGANIZATION 


2  PEO  FIELD  ARTILLERY  SYSTEMS 
SFAE  FAS  PM 

H  GOLDMAN 
T  MCWILLIAMS 
PICATINNY  ARSENAL  NJ 
07806-5000 

6  PM  SADARM 
SFAE  GCSS  SD 
COLB  ELLIS 
M  DEVINE 
R  KOWALSKI 
W  DEMASSI 
J  PRITCHARD 
SHROWNAK 
PICATINNY  ARSENAL  NJ 
07806-5000 

1  COMMANDER 

US  ARMY  ARDEC 
PRODUCTION  BASE 
MODERN  ACTY 
AMSMCPBMK 
PICATINNY  ARSENAL  NJ 
07806-5000 

3  COMMANDER 

U  S  ARMY  TACOM 
PM  TACTICAL  VEHICLES 
SFAETVL 
SFAETVM 
SFAETVH 

6501  ELEVEN  MILE  RD 
WARREN  MI  48397-5000 

1  COMMANDER 

US  ARMY  TACOM 
PM  ABRAMS 
SFAE  ASM  AB 
6501  ELEVEN  MILE  RD 
WARREN  MI  48397-5000 

1  COMMANDER 

U  S  ARMY  TACOM 

PMBFVS 

SFAE  ASM  BV 

6501  ELEVEN  MILE  RD 

WARREN  MI  48397-5000 
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1 

2 


1 

1 

1 


ORGANIZATION 

NO.  OF 
COPIES 

ORGANIZATION 

COMMANDER 

1 

COMMANDER 

USARMYTACOM 

USARMYTACOM 

PM  AFAS 

PM  GROUND  SYSTEMS 

SFAE  ASM  AF 

INTEGRATION 

6501  ELEVEN  MILE  RD 

SFAE  GCSS  W  GSI 

WARREN  MI  48397-5000 

R  LABATILLE 

COMMANDER 

6501  ELEVEN  MILE  RD 
WARREN  MI  48397-5000 

U  S  ARMY  TACOM 

PM  SURV  SYS 

1 

COMMANDER 

SFAE  ASM  SS 

USARMYTACOM 

TDEAN 

CHIEF  ABRAMS  TESTING 

SFAE  GCSS  W  GSI M 

SFAE  GCSS  W  AB  QT 

D COCHRAN 

T  KRASKIEWICZ 

6501  ELEVEN  MILE  RD 

6501  ELEVEN  MILE  RD 

WARREN  MI  48397-5000 

WARREN  MI  48397-5000 

COMMANDER 

1 

COMMANDER 

U  S  ARMY  TACOM 

US  ARMY  TACOM 

PM  RDT&E 

AMSTA  SF 

SFAE  GCSS  WAB 

WARREN  MI  48397-5000 

J  GODELL 

6501  ELEVEN  MILE  RD 

1. 

COMMANDER 

WARREN  MI  48397-5000 

SMCWV  QAE  Q 

B  VANINA 

BLDG  44 

WATER  VLIET  ARSENAL 
WATERVLIET  NY  12189-4050 

14  COMMANDER 

USARMYTACOM 

ASMTATRR 


J  CHAPIN 

COMMANDER 

R  MCCLELLAND 

U  S  ARMY  TACOM 

D  THOMAS 

PMBFV 

J BENNETT 

SFAE  GCSS  WBV 

D  HANSEN 

S  DAVIS 

AMSTA  JSK 

6501  ELEVEN  MILE  RD 

S  GOODMAN 

WARREN  MI  48397-5000 

J  FLORENCE 

K  IYER 

COMMANDER 

J  THOMSON 

USARMYTACOM 

AMSTA  TR  D 

PM  LIGHT  TACTICAL 

D  OSTBERG 

VEHICLES 

L  HINOJOSA 

AMSTA  TR  S 

BRAJU 

AJJ  MILLS  MS  209 

AMSTA  CSSF 

6501  ELEVEN  MILE  RD 

H  HUTCHINSON 

WARREN  MI  48397-5000 

F  SCHWARZ 

WARREN  MI  48397-5000 

COMMANDER 
U  S  ARMY  TACOM 
PM  SURVTVABLE  SYSTEMS 
SFAEGCSSWGSIH 
MRYZYI 

6501  ELEVEN  MILE  RD 
WARREN  MI  48397-5000 
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NO.  OF 

COPIES  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


1  COMMANDER 
SMCWV  SPM 

T  MCCLOSKEY 
BLDG  253 

WATER  VLIET  ARSENAL 
WATERVLIET  NY  12189-4050 

10  BENET  LABORATORIES 
AMSTA  AR  CCB 
R  FISCELLA 
GD  ANDREA 
M  SCAVULO 
G  SPENCER 
P  WHEELER 
K  MINER 
J  VASILAKIS 
G  FRIAR 
R  HASENBEIN 
SMCAR  CCB  R 
S  SOPOK 

WATERVLIET  NY  12189 

2  TSM  ABRAMS 
ATZKTS 
SJABURG 

W  MEINSHAUSEN 
FT  KNOX  KY  40121 

3  ARMOR  SCHOOL 
ATTN  ATZK  TD 
RBAUEN 
JBERG 
APOMEY 

FT  KNOX  KY  40121 

2  HQ  IOC  TANK  AMMO  TEAM 

AMSIO  SMT 
R  CRAWFORD 
W  HARRIS 

ROCK  ISLAND  IL  61299-6000 

1  DIRECTOR 

U  S  ARMY  AMCOM 
SFAE  AV  RAM  TV 
D  CALDWELL 
BUILDING  5300 

REDSTONE  ARSENAL  AL  35898 


4  DIRECTOR 

US  ARMY  CECOM 
NIGHT  VISION  AND  ELECTRONIC 
SENSORS  DIRECTORATE 
AMSEL  RD  NV  CM  CCD 
R  ADAMS 
R  MCLEAN 
AYINGST 
AMSEL  RD  NV  VISP 
E  JACOBS 
10221  BURBECK  RD 
FT  BELVOIR  VA  22060-5806 

2  CDR  USA  AMCOM 

AVIATION  APPLIED  TECH  DIR 
JSCHUCK 

FT  EUSTIS  VA  23604-5577 

1  U  S  ARMY  CRREL 
PDUTTA 
72  LYME  RD 
HANOVER  NH  03755 

1  US  ARMYCERL 
RLAMPO 

2902  NEWMARK  DR 
CHAMPAIGN  IL  61822 

2  US  ARMY  CORP  OF  ENGINEERS 
CERDCTLIU 
CEWETTTAN 

20  MASS  AVE  NW 
WASHINGTON  DC  20314 

10  DIRECTOR 

US  ARMY  NATL  GRND  INTEL  CTR 
D  LETTER 
SETIELMAN 
M  HOLTUS 
M  WOLFE 
S  MINGLEDORF 
H  C  ARDLEIGH 
J  GASTON 

W  GSTATTENBAUER 
R  WARNER 
J  CRIDER 

220  SEVENTH  STREET  NE 
CHARLOTTESVILLE  VA  22091 
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COPIES  ORGANIZATION 


6  US  ARMY  SBCCOM 

SOLDIER  SYSTEMS  CENTER 
BALLISTICS  TEAM 
J  WARD 

MARINE  CORPS  TEAM 

J  MACKIE  WICZ 

BUS  AREA  ADVOCACY  TEAM 

W  HASKELL 

SSCNC  WST 

W  NYKVIST 

TMERREL 

S  BEAUDOIN 

KANSAS  ST 

NATICK  MA  01760-5019 

1  US  ARMY  COLD 

REGIONS  RSCH  & 

ENGRNG  LAB 
PDUTTA 
72  LYME  RD 
HANOVER  NH  03755 

1  SYSTEM  MANAGER  ABRAMS 

BLDG  1002  RM  110 
ATZK  TS  LTC  J  H  NUNN 
FT  KNOX  KY  40121 

9  US  ARMY  RESEARCH  OFFICE 

A  CROWSON 
J CHANDRA 
H EVERETT 
J  PRATER 
R  SINGLETON 
G  ANDERSON 
D  STEPP 
D  KISEROW 
J CHANG 
PO  BOX  12211 

RESEARCH  TRIANGLE  PARK  NC 
27709-2211 

1  DIRECTORATE  OF  CMBT 
DEVELOPMENT 
CKJORO 

320  ENGINEER  LOOP  STE  141 
FT  LEONARD  WOOD  MO  65473-8929 


1  COMMANDANT 

U  S  ARMY  FIELD  ARTILLERY  CTR 
AT  FT  SILL 

ATFS  CD  LTC  BUMGARNER 
FT  SILL  OK  73503  5600 

1  CHIEF  USAIC 

LTC  TJ  CUMMINGS 
ATZB  COM 

FT  BENNING  GA  31905-5800 

1  NAVAL  AIR  SYSTEMS  CMD 
J  THOMPSON 
48142  SHAW  RD  UNIT  5 
PATUXENT  RIVER  MD  20670 

1  NAVAL  SURFACE  WARFARE  CTR 
DAHLGREN  DIV  CODE  G06 
DAHLGRENVA  22448 

1  NAVAL  SURFACE  WARFARE  CTR 
TECH  LIBRARY  CODE  323 
17320  DAHLGREN  RD 
DAHLGRENVA  22448 

3  NAVAL  RESEARCH  LAB 
I WOLOCK  CODE  6383 
R  BADALIANCE  CODE  6304 
LGAUSE 

WASHINGTON  DC  20375 

1  NAVAL  SURFACE  WARFARE  CTR 
CRANE  DIVISION 

M  JOHNSON  CODE  20H4 
LOUISVILLE  KY  40214-5245 

2  COMMANDER 

NAVAL  SURFACE  WARFARE  CTR 
CADEROCK  DIVISION 
R  PETERSON  CODE  2020 
M  CRITCHFIELD  CODE  1730 
BETHESDAMD  20084 

2  NAVAL  SURFACE  WARFARE  CTR 
U  SORATHLA 
C  WILLIAMS  CD  6551 
9500  MACARTHUR  BLVD 
WEST  BETHESDA  MD  20817 
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NO.  OF 
COPIES 

1 

2 

1 

8 


1 

1 

1 

2 


ORGANIZATION 


DAVID  TAYLOR  RESEARCH  CTR 
SHIP  STRUCTURES  &  PROTECTION 
DEPARTMENT  CODE  1702 
J  CORRADO 
BETHESDA  MD  20084 

DAVID  TAYLOR  RESEARCH  CTR 
R  ROCKWELL 
W  PHYILLAIER 
BETHESDA  MD  20054-5000 

OFFICE  OF  NAVAL  RESEARCH 
D  SIEGEL  CODE  351 
800  N  QUINCY  ST 
ARLINGTON  VA  22217-5660 

NAVAL  SURFACE  WARFARE  CTR 
J  FRANCIS  CODE  G30 
D  WILSON  CODE  G32 
RD  COOPER  CODE  G32 
J  FRAYSSE  CODE  G33 
E  ROWE  CODE  G33 
T  DURAN  CODE  G33 
L  DE  SIMONE  CODE  G33 
R  HUBBARD  CODE  G33 
DAHLGREN  VA  22448 

NAVAL  SEA  SYSTEMS  CMD 
DLIESE 

2531  JEFFERSON  DAVIS  HIGHWAY 
ARLINGTON  VA  22242-5160 

NAVAL  SURFACE  WARFARE 
M  LACY  CODE  B02 
17320  DAHLGREN  RD 
DAHLGREN  VA  22448 

OFFICE  OF  NAVAL  RES 
J KELLY 

800  NORTH  QUINCEY  ST 
ARLINGTON  VA  22217-5000 

NAVAL  SURFACE  WARFARE  CTR 
CARDEROCK  DIVISION 
R  CRANE  CODE  2802 
C  WILLIAMS  CODE  6553 
3A  LEGGETT  CIR 
BETHESDA  MD  20054-5000 


NO.  OF 

COPIES  ORGANIZATION 


1  NAVSEA  OJRI 

PEO  DD21  PMS500 
G  CAMPONESCHI 
2351  JEFFERSON  DAVIS  HWY 
ARLINGTON  VA  22242-5165 

1  EXPEDITIONARY  WARFARE 

DIVN85 
FSHOUP 

2000  NAVY  PENTAGON 
WASHINGTON  DC  20350-2000 

1  AFRLMLBC 

2941  P  STREET  RM  136 
WRIGHT  PATTERSON  AFB  OH 
45433-7750 

1  AFRLMLSS 
R THOMSON 

2179  12TH  STREET  RM  122 
WRIGHT  PATTERSON  AFB  OH 
45433-7718 

2  AFRL 

F ABRAMS 
J  BROWN 
BLDG  653 

2977  P  STREET  STE  6 
WRIGHT  PATTERSON  AFB  OH 
45433-7739 

1  AFRL  MLS  OL 

7278  4TH  STREET 
BLDG  100  BAY  D 
L COULTER 

HILL  AFB  UT  84056-5205 
1  OSD 

JOINT  CCD  TEST  FORCE 
OSD  JCCD  R  WILLIAMS 
3909  HALLS  FERRY  RD 
VICKSBURG  MS  29180-6199 

1  DEFENSE  NUCLEAR  AGENCY 

INNOVATIVE  CONCEPTS  DIV 
RROHR 

6801  TELEGRAPH  RD 
ALEXANDRIA  VA  22310-3398 
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NO.  OF 
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1 

3 

2 

1 

2 

1 

1 


ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


WATERWAYS  EXPERIMENT 
D SCOTT 

3909  HALLS  FERRY  RDSCC 
VICKSBURG  MS  39180 

DARPA 

M  VANFOSSEN 
S  WAX 

L  CHRISTODOULOU 
3701  N  FAIRFAX  DR 
ARLINGTON  VA  22203-1714 

SERDP  PROGRAM  OFC  PM  P2 
C  PELLERIN 
B  SMITH 

901  N  STUART  ST  SUITE  303 
ARLINGTON  VA  22203 

FAA 

MIL  HDBK  17  CHAIR 
LILCEWICZ 
1601  LIND  AVE  SW 
ANM  115N 
RENTON  VA  98055 

FAA  TECH  CENTER 
D  OPLINGER  AAR  431 
P  SHYPRYKEVICH  AAR  431 
ATLANTIC  CITY  NJ  08405 

OFC  OF  ENVIRONMENTAL  MGMT 
US  DEPT  OF  ENERGY 
P  RITZCOVAN 
19901  GERMANTOWN  RD 
GERMANTOWN  MD  20874-1928 

LOS  ALAMOS  NATL  LAB 
F  ADDESSIO 
MS  B216 
PO  BOX  1633 
LOS  ALAMOS  NM  87545 

OAK  RIDGE  NATL  LAB 
R  M  DAVIS 
PO  BOX  2008 

OAK  RIDGE  TN  37831-6195 


5  DIRECTOR 

LAWRENCE  LIVERMORE 
NATL  LAB 
R  CHRISTENSEN 
S  DETERESA 
F  MAGNESS 
M  FINGER  MS  313 
M  MURPHY  L  282 
PO  BOX  808 
LIVERMORE  CA  94550 

7  NIST 

RPARNAS 

JDUNKERS 

M  VANLANDINGHAM  MS  8621 
J  CHIN  MS  8621 
D  HUNSTON  MS  8543 
J  MARTIN  MS  8621 
D  DUTHINH  MS  861 1 
100  BUREAU  DR 
GAITHERSBURG  MD  20899 

1  OAK  RIDGE  NATL  LAB 

C  EBERLE  MS  8048 
PO  BOX  2009 
OAK  RIDGE  TN  37831 

1  OAK  RIDGE  NATL  LAB 

CD  WARREN  MS  8039 
PO  BOX  2009 
OAK  RIDGE  TN  37922 

4  DIRECTOR 

SANDIA  NATL  LABS 
APPLIED  MECHANICS  DEPT 
DIVISION  8241 
W  KAWAHARA 
KPERANO 
D  DAWSON 
PNIELAN 
PO  BOX  969 

LIVERMORE  CA  94550-0096 

1  LAWRENCE  LIVERMORE 
NATIONAL  LAB 
M  MURPHY 
PO  BOX  808  L  282 
LIVERMORE  CA  94550 
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3  NASA  LANGLEY  RESEARCH  CTR 

MS  266 
AMSRL  VS 
W  ELBER 
F  BARTLETT  JR 
G FARLEY 

HAMPTON  VA  23681-0001 

1  NASA  LANGLEY  RESEARCH  CTR 

T  GATES  MS  188E 
HAMPTON  VA  23661-3400 

1  USDOT  FEDERAL  RAILROAD 

RDV  31  M  FATEH 
WASHINGTON  DC  20590 

1  DOT  FHWA 

JSCALZI 

400  SEVENTH  ST  SW 
3203  HNG  32 
WASHINGTON  DC  20590 

1  FHWA 

EMUNLEY 

6300  GEORGETOWN  PIKE 
MCLEAN  VA  22101 

1  CENTRAL  INTELLIGENCE  AGENCY 

OTI WDAG  GT 
WL  WALTMAN 
PO  BOX  1925 
WASHINGTON  DC  20505 

1  MARINE  CORPS  INTEL  ACTY 

D  KOSITZKE 

3300  RUSSELL  RD  SUITE  250 
QUANnCO  VA  22134-5011 

1  NATL  GRND  INTELLIGENCE  CTR 

DIRECTOR 
IANGTMT 
220  SEVENTH  ST  NE 
CHARLOTTESVILLE  VA  22902-5396 

1  DIRECTOR 

DEFENSE  INTELLIGENCE  AGENCY 
TA  5  K  CRELLING 
WASHINGTON  DC  20310 


NO.  OF 

COPIES  ORGANIZAnON 


1  GRAPHITE  MASTERS  INC 

J  WILLIS 

3815  MEDFORD  ST 

LOS  ANGELES  CA  90063-1900 

1  ADVANCED  GLASS  FIBER  YARNS 

T  COLLINS 

281  SPRING  RUN  LN  STE  A 
DOWNINGTON  PA  19335 

1  COMPOSITE  MATERIALS  INC 

D  SHORTT 
19105  63  AVENE 
PO  BOX  25 

ARLINGTON  WA  98223 

1  COMPOSITE  MATERIALS  INC 

R  HOLLAND 
11  JEWEL  COURT 
ORINDA  CA  94563 

1  COMPOSITE  MATERIALS  INC 
C  RILEY 

14530  S  ANSON  AVE 
SANTA  FE  SPRINGS  CA  90670 

2  COMPOSIX 
D BLAKE 

L  DIXON 

120  0  NEILL  DR 

HEBRUN  OHIO  43025 

4  CYTEC  FEBERITE 
R DUNNE 
DKOHLI 
MGILLIO 
RMAYHEW 
1300  REVOLimON  ST 
HAVRE  DE  GRACE  MD  21078 

2  SIMULA 

J  COLTMAN 
RHUYETT 
10016  S51ST  ST 
PHOENIX  AZ  85044 

1  SIOUX  MFG 
BKRIEL 
PO  BOX  400 
FT  TOTTEN  ND  58335 
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COPIES  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


2  PROTECTION  MATERIALS  INC 
M  MILLER 

F  CRILLEY 
14000  NW  58  CT 
MIAMI  LAKES  FL  33014 

3  FOSTER  MILLER 
J  J  GASSNER 

M  ROYLANCE 
WZUKAS 
195  BEAR  HILL  RD 
WALTHAM  MA  02354-1 196 

1  ROM  DEVELOPMENT  CORP 
ROMEARA 

136  SWINEBURNE  ROW 
BRICK  MARKET  PLACE 
NEWPORT  RI 02840 

2  TEXTRON  SYSTEMS 
T FOLTZ 

M TREASURE 
201  LOWELL  ST 
WILMINGTON  MA  08870-2941 

1  JPS  GLASS 

L CARTER 
PO  BOX  260 
SLATER  RD 
SLATER  SC  29683 

1  O  GARA  HESS  &  EISENHARDT 
M  GILLESPIE 

9113  LESAINT  DR 
FAIRFIELD  OH  45014 

2  MILLIKEN  RESEARCH  CORP 
HKUHN 

M  MACLEOD 
PO  BOX  1926 
SPARTANBURG  SC  29303 

1  CONNEAUGHT  INDUSTRIES  INC 

J SANTOS 
PO  BOX  1425 
COVENTRY  RI  02816 


1  BATTELLE 

CR  HARGREAVES 
505  KING  AVE 
COLUMBUS  OH  43201-2681 

2  BATTELLE  NATICK  OPERATIONS 
J  CONNORS 

BHALPIN 

209  W  CENTRAL  ST 

STE  302 

NATICK  MA  01760 

1  BATTELLE  NW  DOE  PNNL 

T  HALL  MS  K231 
BATTELLE  BLVD 
RICHLAND  WA  99352 

3  PACIFIC  NORTHWEST  LAB 
M  SMITH 

G  VAN  ARSDALE 
R  SHIPPELL 
PO  BOX  999 
RICHLAND  WA  99352 

1  ARMTEC  DEFENSE  PRODUCTS 
SDYER 

85  901  AVE  53 
PO  BOX  848 
COACHELLA  CA  92236 

2  ADVANCED  COMPOSITE 
MATLS  CORP 
PHOOD 

J  RHODES 

1525  S  BUNCOMBE  RD 
GREER  SC  29651-9208 

2  GLCCINC 
JRAY 

M  BRADLEY 

103  TRADE  ZONE  DR 

STE26C 

WEST  COLUMBIA  SC  29170 

2  AMOCO  PERFORMANCE 
PRODUCTS 
M  MICHNO  JR 
J  BANISAUKAS 
4500  MCGINNIS  FERRY  RD 
ALPHARETTA  GA  30202-3944 
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COPIES  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


1  SAIC 

M  PALMER 

2109  AIR  PARK  RD  S  E 

ALBUQUERQUE  NM  87106 

1  SAIC 

ATTN  G  CHRYSSOMALLIS 
3800  W  80TH  ST  STE  1090 
BLOOMINGTON  MN  55431 

1  AAI  CORPORATION 

DR  T  G  STASTNY 
PO  BOX  126 

HUNT  VALLEY  MD  21030-0126 

1  JOHN  HEBERT 

PO  BOX  1072 

HUNT  VALLEY  MD  21030-0126 

12  ALLIANT  TECHS  YSTEMS  INC 
CCANDLAND 
C  AAKHUS 
R  BECKER 
B  SEE 

N VLAHAKUS 
RDOHRN 
SHAGLUND 
D  FISHER 
W  WORRELL 
R  COPENHAFER 
M  HISSONG 
DKAMDAR 
600  2ND  ST  NE 
HOPKINS  MN  55343-8367 

3  ALLIANT  TECHS  YSTEMS  INC 

J  CONDON 
ELYNAM 
J GERHARD 
WV01  16  STATE  RT  956 
PO  BOX  210 

ROCKET  CENTER  WV  26726-0210 

1  APPLIED  COMPOSITES 

WGRISCH 

333  NORTH  SIXTH  ST 
ST  CHARLES  IL  60174 


1  PROJECTILE  TECHNOLOGY  INC 

515  GILES  ST 

HAVRE  DE  GRACE  MD  21078 

1  CUSTOM  ANALYTICAL 
ENG  SYS  INC 

A  ALEXANDER 
13000  TENSOR  LN  NE 
FLINTSTONE  MD  21530 

2  LORAL  VOUGHT  SYSTEMS 
G  JACKSON 

KCOOK 

1701  W  MARSHALL  DR 
GRAND  PRAIRIE  TX  7505 1 

5  AEROJET  GEN  CORP 

D  PILLASCH 
T COULTER 
C  FLYNN 
DRUBAREZUL 
M  GREINER 

1 100  WEST  HOLLYVALE  ST 
AZUSA  CA  91702-0296 

3  HEXCELINC 
RBOE 

F  POLICELLI 
JPOESCH 
PO  BOX  98 
MAGNA  UT  84044 

3  HERCULES  INC 

G  KUEBELER 
J  VERMEYCHUK 
B  MANDERVILLE  JR 
HERCULES  PLAZA 
WILMINGTON  DE  19894 

1  BRIGS  COMPANY 

J  BACKOFEN 
2668  PETERBOROUGH  ST 
HERDON  VA  22071-2443 

1  ZERNOW  TECHNICAL  SERVICES 

LZERNOW 

425  W  BONITA  AVE  STE  208 
SAN  DIMAS  CA  91773 
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COPIES  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


2  OLIN  CORPORATION 

FLINCHBAUGHDIV 
E  STEINER 
B  STEWART 
PO  BOX  127 
RED  LION  PA  17356 

1  OLIN  CORPORATION 

L  WHITMORE 
10101  9TH  ST  NORTH 
ST  PETERSBURG  FL  33702 

1  DOW  UT 

S  TIDRICK 
15  STERLING  DR 
WALLINGFORD  CT  06492 

5  SIKORSKY  AIRCRAFT 
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